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ABSTRACT 

We observed the new X-ray transient and black-hole candidate XTE J1652— 453 si- 
multaneously with XMM-Newton and the Rossi X-ray Timing Explorer (RXTE). The 
C| | observation was done during the decay of the 2009 outburst, when XTE J1652— 453 

was in the hard-intermediate state. The spectrum shows a strong and broad Fe emis- 
sion line with an equivalent width of ~ 450 eV. The profile is consistent with that 
of a line being produced by reflection off the accretion disc, broadened by relativistic 
effects close to the black hole. The best-fitting inner radius of the accretion disc is 
' ~ 4 gravitational radii. Assuming that the accretion disc is truncated at the radius 

of the innermost stable circular orbit, the black hole in XTE J1652— 453 has a spin 
parameter of ~ 0.5. The power spectrum of the RXTE observation has an additional 
variability component above 50 Hz, which is typical for the hard-intermediate state. 
No coherent quasi-periodic oscillations at low frequency are apparent in the power 
spectrum, which may imply that we view the system at a rather low inclination angle. 
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1 INTRODUCTION l200tj ; iBellonil [2010| . and references therein). In the bright 

phase of an outburst, at X-ray luminosities Lx ~ 10 37 -10 39 
• i-H Most binaries with a stellar-mass black- hole primary are X- _i , . , , , , 

k j. i*i el S s > the X-ray spectrum is dominated by thermal 

ray transients m which matter is accreted from the com- . . , , _ , . . r ,. ,,. , 

r > J emission below ~ 5 kev originating from an optically thick 
i pamon star onto the black hole BH . These systems are , . . .. . ,. ,. ... , . . , 

rrt ■ ; / and geometrically thm accretion disc, with a typical tem- 

characterised by a sudden and rapid rise m intensity (reach- r , . T . , J . , , .. 

,,„,,,, perature 01 ~ 1 kev, that extends down to the marginally 
mg a maximum m about one to several weeks followed by . , , ,., nr, — , „ — n Ln-j, ™ . , 

b . a 1 stable orbit (Shakura & Sunyacv 1973). This phase is 

a gradual decline, the so-called outbursts typically last sev- ,, r , . , . , , e . , , , 

generally referred to as the soft or high/sort state, and 

eral months; eventually the source becomes too weak to be . „ , , , , . . . . , .,, 

is further characterised by a low level of variabifity, with 
detected with current X-ray instruments (for typical light . , r , . , , r r . „„ ,, 

„ r- — . . \ , „. ' tX I 1. ° typical fractional rms vafues of a few per cent. When the 

curves, see Fig. 4-9 in lRemillard fe McClinto ck 2006). The , . , , , , , , \ t , ^ 

— — ' luminosity decreases to about 1-4 per cent of the Jiddington 

large excursions m the intensity seen m X-ray transients , . rrz — ni7^7T^l\ ^ i 

luminosity (Maccaronc 2003), the source makes a transition 
are likefy due to changes m the accretion rate onto the BH, . . ,111 1 , n , , ■ , • , al v 

111,. ?.,.., , , , to the so-called hard or low/hard state, m which the X-ray 

triggered by disc instabilities due to hydrogen lomsation (see . . . . , , , , ... . . 

„ rr -11 tj rr— "■ — in— —I, spectrum is dominated by hard power-law like emission 

Section 2.2.1 in lDone et al.1120071 ; lLasotall200ll ). , . , . y , . 1 . , , , _ , _ „, 

— : — — — : — ' with a typical power-law photon index of 1.5-1.7. the 

During an outburst, not only the intensity, but , . , ' . . . , , , „ 

J high-energy emission is produced by inverse Compton 

also the X-ray spectrum and X-ray variability changes ... „ , , , ,,. „ , , , . 

r— ; — : ' •' , r T „ 1 1 -1 1 , , , * „ " , 1 scattering of photons by a popufation of hot cfectrons m 

jMendez fe van der Khsl 119971 : [McChntock & Rcinillard . , , . F7Z 1 p ■ J ITT^K 

' a region referred to as the corona |lhornc & Price 1975). 

The hard state is further characterised by a high level of 

i, „ ., , . . . , variability, with rms amplitudes of > 20 per cent, with 

lj-mail: niemstraSSastro. rug.nl ■" 1 1- ; 
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strong quasi-perio dic oscillations (QPO) a t low frequencies 
(~ 0.1-10 Hz; e.g., iHoman fc Bellonill2005l ). 

X-ray spectra of BH transients sometimes show broad 
Fe em ission lines at around 6-7 keV (for a review see iMiller! 

2007) . These lines, from fluorescent-Ka emission of neutral 
or ionised Fe, are thought to be produced through irradia- 
tion of the disc by the source of hard photons. The Fe line 
is b roadened by relativistic effects and becomes asymmet- 
ric l|Fabian et al.lll989l ; lLaorH l991). and its profile provides 
information on the disc properties, such as the inclination 
angle, the emissivity, and the inner ra dius. The latter can 
be us ed to measure the BH spin (e.g., iRevnolds fc Fabianl 

2008) , assuming that the disc extends down to the innermost 
stable circular orbit (ISCO) in the soft states. 

In this paper we report on the timing and spectral 
analysis of a contemporaneous XMM-Newton and RXTE 
observation of the new X-ray transient XTE J1652— 453. 
This is the first (and so far the only) time that this source 
has been observed with an instrument with a good spec- 
tral resolution to study the Fe line which was previously 
detected in RXTE Propo rtianal Counter Array (PCA) data 
l|Markwardt et al. 20 09a). Furthermore, the broadband cov- 
erage of the simultaneous XMM-Newton and RXTE data al- 
lows to properly model the underlying continuum spectrum. 
In Section [2] we describe the outburst evolution from RXTE 
observations, we give details about the XMM-Newton and 
RXTE observations, and we describe the steps that we fol- 
lowed to generate the power and energy spectra. We continue 
with the data analysis in Section [3] where we test different 
line and continuum models, and where we consider several 
broadening mechanisms. In this section we also present our 
results from the spectral and timing analysis separately. In 
Section U we discuss our results in the context of the spec- 
tral state, disc geometry, and inferred black hole spin. In the 
appendix we discuss how the choice of background spectrum 
and abundances influences the continuum and line parame- 
ters inferred from the spectrum. 

1.1 The 2009 outburst of XTE J1652 453 

XTE J1652-453 was discovered i n the R XTE /PCA Galac - 
tic bulge scan on June 28 2009 l|Markwardt et al. 2009a! '). 
Follow-up observations with the Swift X-ray Telescope 
(XRT) and RXTE /PCA measured a spectrum consistent 
with a stellar-mass black hole in the canonical high/soft 
state (|Markwardt fc Swank 20091 ; iMarkwardt et al. 2009bh . 
with the spectrum below 10 keV dominated by emission from 
an accretion disc with a temperature of ~ 0.6 keV. The 
RXTE /PCA spectrum also showed evidence of a relatively 
strong emission line at 6.4 keV, most likely originating from 
an Fe-Ka line. However, the line could not be studied in 
detail because of the limited spectral resolution of RXTE. 

A near infrared (NIR) imaging observation of the field 
of XTE J1652-453 on July 11 200 9, provided a possibl e 
NIR counterpart of XTE J1652-453 l|Revnolds et al.ll2009h . 
However, from NIR observations on July 15 and August 
28 2009, the suggested NIR counterpart showed no signifi- 
cant variability, suggesting that the prev ious detection migh t 
have been an unrelated interloper star (|Torres et ali i2009). 
Radio observations with the Australia Telescope Compact 
Array on July 6, 7, and 13 2009, showed emission that may 
be indicative of the decay of an optically thin synchrotron 



event associated w ith the activation of XTE J1652— 453 
l|Calvelo et alj|2009h . 

On August 22 2009, we triggered an XMM-Newton ob- 
servation of XTE J1652-453, simultaneously with RXTE, 
to study the Fe emission line and the disc properties in the 
decline of the outburst. In this paper we present the results 
of these observations. 

Around September 10 2009, the Burst Alert Telescope 
(BAT) on-board Swift showed an increase (by a factor 
of 1.5-2) in the hard X-ray count rate (15-50 keV) of 
XTE J1652— 453, in dicating a transition to the hard state 
i|Coriat et al.l 120091) . On September 14 2009, the RXTE 
spectrum was still consistent with a disc blackbody (kT ~ 
0.6 keV), a power law (index of ~ 1.8), and a line feature at 
~ 6 keV. Later RXTE observations did not r equire a disc 
comp onent, and no line feature was apparent (|Coriat et al.l 
l2009f ). 



2 OBSERVATIONS & DATA REDUCTION 

2.1 Hardness and intensity as function of time 

Since its discovery, XTE J1652— 453 was observed with 
RXTE every - 3-4 days. There are in total 55 RXTE obser- 
vations of XTE J1652— 453, each of them lasting from ~ 1 
to several ks. 

We used the 16-s time-resolution St andard 2 data from 
the PCA (for instrument information see lJahoda et al ]|2006l ) 
to calculate the hard colour and intensity of all 55 observa- 
tions. For each of the five PCA units (PCUs), and for each 
time interval of 16 s, we calculated the intensity, defined as 
the count rate in the 2-20 keV band, and the hard colour, 
defined as the count rate in the 16-20 keV range divided by 
the count rate in the 2-6 keV range. 

To correct for differences in effective area between the 
PCUs and for time-dependent ga in changes of the de t ectors , 
we used the method described in Ivan Straaten et"al] l|2003l ). 
As we did for XTE J1652-453, for each PCU we calculated 
the hard colour and intensity of the Crab nebula, which 
are assumed to be constant, and averaged the 16-s Crab 
colours and intensities per day. For each PCU, we divided 
the 16-s hard colours and intensities of XTE J1652— 453 by 
the corresponding averaged Crab values that were closest in 
time. Finally, we averaged the 16-s colours and intensities 
over all PCUs per observation. 

We further used the 15-50 keV fluxes of XTE 
J1652— 453 measured by Swift /BAT during the outburst. 
These data are provided by the Swift /BAT team and 
are publicly available on the Swift /BAT transient monitor 
pagefl 

2.2 Energy spectrum of XTE J1652-453 

XTE J1652-453 was simultaneously observed with XMM- 
Newton and RXTE on August 22 2009; both observations 
started at the same time (03:39 UTC). The XMM-Newton 
observation (obsID 0610000701) lasted for ~ 42 ks, while 
the RXTE observation (obsID 94432-01-04-00) lasted for 
~ 14 ks. To achieve a broadband spectral coverage, we 

1 http : //swift . gsf c .nasa. gov/docs/swift /results/transients/ 
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extracted the energy spectra of XTE J1652— 453 from the 
XMM-Newton data as well as from the RXTE observation. 
We describe the data reduction of the individual instruments 
in the following sections. 

2.2.1 XMM-Newton 

For the data reduction of the XMM-Newton Observational 
Data Files (ODF) we used the Science Analysis Software 
(SAS) version 9.0, applying the latest calibration files. We 
generated concatenated and calibrated event lists for the 
different instruments, using the standard SAS tasks. The 
European Photon Imaging Camera (EPIC) PN and MOS2 
were operated in timing mode, while MOS1 was operated in 
imaging mode. The Reflection Grating Spectrometers (RGS) 
were set in normal spectroscopy mode. 

For the EPIC data we applied the standard filtering 
criteria, as described in the science threads provided by the 
XMM-Newton team, and we additionally filtered out the 
flaring particle background, leading to a net exposure of 
~ 35 ks. Except for the flaring background, the light curve 
was rather constant during the whole observation. Since 
MOS1 was operated in imaging mode, it was heavily af- 
fected by photon pile-up. After excising the pile-up affected 
central region of the source, the statistics of the MOS1 spec- 
trum were strongly reduced, and we therefore did not include 
this spectrum in the rest of our analysis. Using the task 
epatplot, we evaluated the pile-up fraction in MOS2 and 
PN data, and we found that no photon pile-up was appar- 
ent for MOS2 and PN. To further che ck whether the spec- 
trum was affected by pil e-up (see e.g.. lYamada et ah 2009; 
iDone fc Dfaz Tried 1 2009 ) , we extracted the source spectra 
excluding the central three columns of the point spread func- 
tion (PSF). There are no significant differences between the 
spectrum that includes the core of the PSF and the one that 
excludes it. 

Following the standard procedure, we extracted source 
and background spectra for MOS2 and PN data, and we cre- 
ated redistribution matrices and ancillary files using the SAS 
tools RMFGEN and ARFGEN, respectively. The MOS2 source 
spectrum was extracted using the RAWX columns in 288- 
328, and for single events only (pattern 0). We extracted 
the MOS2 background from a source-free circular region on 
the outer CCDs, ~ 15 arcmin away from the source, using 
pattern 0,1, and 3. 

The standard correction for charge-transfer inefficiency 
(CTI) does not work properly for PN operated in fast mode 
(burst or timing) and bright sources, which ca n lead to a 
gain shift and may affect the energy spectrum (Sal a et all 
I2006T ). Therefore, before extracting the PN spectra, we used 
the new SAS task epfast to correct for possible CTI effects. 
(More information about CTI correction can be found at: 
http : //xmm. esac . esa. int/external/xmm_calibration/ . ) 
We extracted the PN source spectrum using RAWX in 
[30:46]; for the PN background we initially used the RAWX 
columns in [5:21]. However, the background spectrum 
appeared to be contaminated with source photons (see 
Fig. IA7|) . We tried different background extraction regions, 
but since the wings of the PSF extends beyond the CCD 
boundaries, the background spectra taken from the same 
CCD were always contaminated with source photons. 
Using this contaminated background would lead to an 



over-subtraction of the source spectrum, and may affect the 
results of the spectral analysis. Since during our observation 
the source count rate was significantly higher than that of 
the background (source ~ 100 cts s" 1 ; no contaminated 
background < 0.5 cts s -1 ), we decided not to s ubtract any 
background (see also IDone k, Dfaz Trigo 2009). However, 
we investigated whether the best-fitting parameters are 
affected by applying different background corrections. We 
therefore performed a number of spectral fits using different 
background models. Details about these fits can be found 
in the Appendix. 

The EPIC spectra obtained following the standard pro- 
cedures oversample the intrinsic resolution of the detectors 
by a factor of ~ 10. We therefore rebinned the EPIC spec- 
tra with the tool PHARBN0, to have 3 energy channels per 
resolution element, and at least 15 counts per channel. We 
did not apply any systematic error to the EPIC spectra. 

We reduced the RGS data and produced all necessary 
files following the standard procedure. We included the RGS 
spectra at low energies (below 2 keV), only using the first 
order spectra of RGS1 and RGS2. Unless mentioned differ- 
ently, we rebinned the energy channels above ~ 0.9 keV by 
a factor of 10, and the channels below that energy such that 
there are at least 15 counts per bin. No systematic error was 
applied to the RGS data. 



2.2.2 RXTE 

For the RXTE data reduction we used the heasoft tools 
version 6.7. We followed the standard procedure suggested 
by the RXTE tearajj. After standard screening and filter- 
ing we obtained a net exposure of ~ 13 and ~ 4 ks for 
PCA and High Energy X-ray Timing Experiment (HEXTE), 
respectively. Using the tool SAEXTRCT, we extracted the 
PCA spectrum from Standard-2 data, using only events from 
PCU2, being this the best-calibrated detector. We extracted 
the HEXTE spectra using cluster-B events only. The PCA 
background was estimated and the HEXTE background was 
measured with the tools PCABACKEST and hxtback, respec- 
tively. We generated PCA and HEXTE instrument response 
files, using PCARSP and hxtrsp, respectively. Within the 
used energy range the energy channels of the RXTE data 
contained more than 15 counts per bin, such that \ 2 statis- 
tics can be applied, and no channel rebinning was required. 
No systematic errors were added to the HEXTE data, while 
we used 0.6 per cent systematic errors for the PCA data. 



2.3 Power spectrum of XTE J1652 453 

We used the RXTE /PCA data to calculate power density 
spectra of all available observations of XTE J1652— 453. For 
this we used the Event-mode data with 125- /us time res- 
olution over the nominal 2-60 keV range, without any en- 
ergy selection, using data from all available PCUs. For all 55 
RXTE pointed observations, we first produced light curves 
at a 1-s time resolution to identify and exclude drop outs in 
the data. We then divided each observation into segments 



2 http : //virgo .bitp .kiev.ua/docs/xmm_sas/Pawel/reduction/pharbn 

3 http : //heasarc . gsf c .nasa. gov/docs/xte/xtegof .html 
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of continuous 128 s, and we calculated the Fourier trans- 
form of each segment preserving the available time resolu- 
tion. Therefore, each power spectrum extends from 1/128 to 
4096 Hz, with a frequency resolution of 1/128 Hz. Finally, 
we averaged all the 128-s power spectra within an obser- 
vation to produce an averaged power density spectrum per 
observation. 

For each observation we subtracted the de ad-time mod- 
ified P oisson level following the prescription in lZhang et al.l 
|l995l ). Using the corresponding source and background 
count rates, we re-normalised the averaged power sp ectra 
of each observation to rms units l|van der Klis 1989]). We 
then calculated the integrated fractional rms amplitude for 
each observation in the frequency range 1/128-100 Hz. 

In addition, for the RXTE observation (ObsID 94432- 
01-04-00) that was coincident with the XMM-Newton obser- 
vation, we fit the power spectrum with a function consistin g 
of the sum of three Lorentzians (e.g., Belloni et al.l I2OO2T ) . 
The results of this fit are described in Section [3~5l 



3 DATA ANALYSIS AND RESULTS 
3.1 Outburst evolution 

In Fig. [TJ we show the intensity and hard colour as a func- 
tion of time for the outburst of XTE J1652-453. We also 
show the BAT flux and the 1/128-100 Hz fractional rms 
amplit ude. The first detection of this source was on MJD 
55010 (|Markwardt et al. 2009ah . The PCA pointed observa- 
tions started ~ 4 days later, when XTE J1652— 453 had an 
intensity of about 80 mCrab. In the following 2 days the in- 
tensity increased, reaching a maximum of about 140 mCrab 
at MJD ~ 55016. From that point onwards, the intensity 
of the source decayed more or less exponentially until MJD 
~ 55065. After that day the intensity remained approxi- 
mately constant (~ 20 mCrab) for about 55 days, showing 
a dip at around MJD 55080. Starting on MJD 55120, the 
intensity started to decrease again, and ~ 25 days later it 
dropped down below the detection limit of RXTE. In total, 
the outburst lasted approximately 145 days. 

From the beginning of the outburst the hard colour 
increased slowly, roughly linearly with time until MJD ~ 
55065, and at that point it started to increase more rapidly. 
The faster increase occurred at the same time as the inten- 
sity stopped its exponential decay (see above). After this, 
the hard colour continued increasing, showing a small dip 
at MJD 55080, and it further increased by a factor of ~ 3 
within ~ 6 days. The time of the sudden increase in the 
hard colour coincides with the dip in the intensity curve. 
Starting on MJD 55087, and for about 13 days, the hard 
colour remained approximately constant. After MJD 55120 
the hard colour started to increase again, at the time the 
intensity started to decrease again, until XTE J1652— 453 
became too weak to be detected with RXTE. The sudden 
increase in the hard colour at around MJD 55080 also coin- 
cides with a slight increase of the BAT flux light curve, as 
shown in the third panel of Fig. fJJ 

The moment of our XMM-Newton observation is 
marked by a filled black circle in Fig. [T] Our observation 
samples the outburst at the moment that the intensity sta- 
bilised at about ~ 20 mCrab, a day before the first increase 
in the hard colour. 
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Figure 1. From top to bottom: Intensity (2-20 keV), Hard colour 
(16-20 keV/ 2-6 keV), BAT flux (15-50 keV), and fractional rms 
amplitude (1/128-100 Hz, full PCA band), respectively. Intensi- 
ties, hard colours and fractional rms amplitudes are the average 
per day of each pointed PCA observation. The BAT measure- 
ments were taken from the BAT transient monitoring. The black 
filled circle marks the moment of the XMM-Newton observation. 



In the bottom panel of Fig. [T]we show the 1/128-100 
Hz fractional rms amplitude as a function of time during the 
outburst of XTE J1652— 453. In the first few observations 
the rms amplitude is low, at around 2-4 per cent. As the 
outburst progresses, and the source intensity decreases, the 
rms amplitude increases more or less linearly with time, up 
to ~ 30 per cent on MJD 55094, then it decreases slightly 
to ~15-20 per cent for about a week, and then increases 
again up to 30-40 per cent at the end of the outburst. At 
the time of the simultaneous XMM-Newton observation, on 
MJD 55065, the 1/128-100 Hz fractional rms amplitude of 
XTE J1652— 453 was somewhat intermediate, 13.3±1.2 per 
cent. 



3.2 Spectral state at the moment of the 
XMM-Newton observation 

From Fig. [1] it is apparent that at the beginning of the 
outburst XTE J1652— 453 was in the high/soft state, with 
high intensity, low hard colo ur, and low 1/128-100 Hz frac- 
tional rms amplitude (e.g., Me ndez fc van der Klis! Il997i : 
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Figure 2. Hardness intensity diagram of the 2009 outburst of 
XTE J1652-453 (black), and of the 2002/2003 outburst of GX 
339—4 (grey). Intensities and hard colours are given in Crab units 
(for more details see Section [2j) • For the HID of GX 339-4 we 
mark ed the different states identified as in iBelloni et al.l ll2005l l 
and lHoman fe Bellonil l|2005h . We also marked the moment of the 
XMM-Newton observation. 



iHoman fc Bellonill2005l ). Around MJD 55080, there is a clear 
transition to the low/hard state, with a sharp increase of 
the hard colour, while both the BAT flux and the 1/128- 
100 Hz fractional rms amplitude increase slightly, but sig- 
nificantly. Transitions from the high/soft state to the inter- 
medi a te state (soft- or hard-intermediate; IHoman fc Bellonil 
120051 ; iBelloml 120101 ) are generally smooth transitions, and 
without multi-wavelength observations in combination with 
detailed timing analysis it is difficult to say when this tran- 
sition exactly happened. 

In Fig.[2]we show the hardness intensity diagram (HID) 
for the outburst of XTE J1652— 453. For comparison we also 
show the HID of the black-hole candidate GX 339—4 of 
the 2002/2003 outburst. Both diagrams show the intensi- 
ties and hard colours in Crab units, both calculated in the 
same manner as described in Section 12.11 In the HID of 
G X 339—4 we m a rked the di fferent spectral sta t es ide ntified 
bv lBelloni etafl (|2005l ) and lHoman fc Bellonil (|2005l ). Sug- 
gested by Fig. [2] together with the 1/128-100 Hz integrated 
rms amplitude of ~ 15 per cent, which is too high for the 
soft or the soft-intermediate state, we conclude that at the 
moment of the XMM-Newton observation, XTE J1652-453 
was in the hard-intermediate state. 



3.3 Spectroscopy 

For the spectral analysis of XT E J 1652— 453 we used xspec 
version 12.5.1 l|Arnaudl [ 1996). Since XMM-Newton and 
RXTE started their observations at the same time, to- 
gether with the rather constant light curve during the whole 
XMM-Newton observation, we fit the XMM-Newton and 
RXTE spectra simultaneously. We initially used a contin- 
uum model including a power law (powerlaw in xspec) 
and a disc blackbody (diskbb; [Mitsuda et al.lll984l ). We in- 
cluded a multiplicative factor in our fit to account for dif- 



ferences in the absolute flux calibration between the differ- 
ent instruments, and to account for interstellar absorption 
we used the photoelect ric absorption comp onent phabs, as- 
su ming cross sections o f lVerner et alj (Il996t ) and abundances 
of I Wilms et~ai] (|2000T l. (In the append ix we show how the 
choice of cross sections and abundances may affect the pa- 
rameters.) 

A strong and broad emission line at around 6-7 keV 
was apparent in the spectra, most likely due to Fe fluores- 
cence lines (of neutral or ionised Fe), and we initially fit the 
line with a Gaussian. However, a Gaussian line did not pro- 
vide a good fit (xt — 2.6); the Gaussian fit a large part of 
the line between 6-7 keV, but it was not able to fit the red 
wing extending down to 4-5 keV. Guided by the asymmet- 
ric profile, we fit the line with a rela tivistic line profile (e.g., 
for a Kerr- metric: Laor; lLaorlll99ll ). We tried different line 
and continuum components to test whether the line param- 
eters depend on the choice of line and/or continuum model. 
Details and results of this can be found in Section \3. 3. 2 1 

When fitting the XMM-Newton/ RXTE spectra, several 
features appeared in the fitting residuals. Some of them are 
due to calibration issues of XMM-Newton data in timing 
mode, some of them are caused by cross-calibration issues 
between the different instruments, and others might be real. 
We discuss these features individually in the following Sec- 
tion. 

Unless stated differently, all uncertainties given in this 
paper are at a 90 per cent confidence level (A^ 2 = 2.71), and 
upper limits are given at a 95 per cent confidence level. The 
equivalent widths (EW) of emission and absorption features 
are both quoted with positive values. 



3.3.1 Instrumental/ Cross-calibration issues 

Here we list all features and issues that arose during 
our spectral analysis which are most likely related with 
uncertainties in the (cross) calibration of the different 
instruments. 

Excess at ~ 1 keV: 

When we fit the PN spectrum in the 0.7-11 keV range 
with a line and continuum model, phabs(powerlaw + 
diskbb + Laor), we found a strong and broad emission 
feature at ~ 1 keV, extending up to ~ 1.3 keV (see Fig. [3]). 
The ~ 1 keV excess was also apparent in the MOS2 data, 
although the profile of the excess was different. Since both 
detectors showed this excess we cannot discard that it may 
be real, and we considered the possibility tha t it might be 
due t o blurred Fe-L or Ne-X emission lines (see lFabian et all 
2009). We tried to fit this excess with either a Gaussian com- 
ponent or a relativistic line profile, but none of these models 
gave a good fit. No excess at ~ 1 keV was apparent in the 
RGS data, although the lower effective area does not allow 
us to discard it completely. 

Similar excesses at ~ 1 keV have been seen in other 
sources observed with PN operated in timing mode (e.g ., 
iMartocchia et al.ll2006l ; ISala et al.ll2008l ; iBoirin et ai]|2005r i . 
The excess appears in that part of the spectrum where the 
flux drops abruptly due to interstellar absorption, and the 
EW of the excess a ppears to increase as a fu nction of this 
absorption (see e.g., iDone fc Diaz Tried 120091 ). Given this, 
a calibration issue related with the redistribution matrix of 
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Figure 3. An excess at ~ 1 keV is apparent in the PN spectrum 
of XTE J1652— 453. We show the data/model ratio when we fit 
the data with a phabs(powerlaw + diskbb + Laor) model. 
The spectrum is not corrected for background, and instrumental 
residuals are still present. 




Energy (keV) 



Figure 4. Zoom in of the PN spectrum of XTE J1652-453. 
The absorption feature at ~ 9.2 keV is apparent in the lower 
panel, where we show the residuals when fitting the data with the 
PHABS(POWERLAW + DISKBB + Laor) model. In the sub-panel we 
show a zoom in of the residuals to show the apparent absorption 
features at ~ 8.5—9 keV. 



the timing-mode data cannot be excluded at the moment. 
Due to the uncertain origin of this excess (a thorough 
investigation is beyond the scope of this paper), we ignored 
the EPIC data below 1.3 keV and included the RGS data 
to cover the spectrum below that energy. 

Energy -dependent discrepancy between PCA and PN: 

Fitting the PCA data in the 3-20 keV range shows 
an energy-dependent deviation between the PN and PCA 
spectra below ~ 7 keV. (With respect to the PN data, the 
PCA data show an excess that increases with decreasing 
energy.) Fitting the PN and PCA spectra individually (in 
the 1.3-11 keV and 3-20 keV ranges, respectively), there 
is a discrepancy in the power-law photon index, which is 
most likely due to a mismatch in cross-calibration between 
the PN and PCA instruments. To avoid the effects of this 
apparent energy-dependent discrepancy between PN and 
PCA data, in the rest of our analysis we only used PCA 
data in the 7-20 keV range. 

Absorption features: 

The fit residuals of the phabs(powerlaw + diskbb 
+ Laor) model to the EPIC data show two absorption 
features at around 1.8 and 2.1 keV, consistent with the 
Si-K (1.84 keV) and Au-M (2.21 keV) edges, which are 
most likely of instrumental origin. To reduce the impact 
of these lines, we could fit each of them with a Gaus- 
sian absorption component at 1.85 ± 0.01 and 2.13 ± 0.02 
keV, and widths of 55 ± 5 eV and 110 ± 10 eV, respec- 
tively, although some residuals between ~ 1.8-2.2 keV still 
remained, probably due to problems with the gain and 
the CTI correction. In order to keep the model as simple 
as possible, instead of fitting these two instrumental fea- 
tures with absorption components, we ignored all energy 
channels between 1.75-2.25 keV so that these instrumental 
edges will not affect the results of the analysis. However, 
we caution that the edges, if caused by residual problems 
of the CTI correction, indicate that a slight energy shift 



may appear in the whole energy band (see XMM-Newton 
Calibration Documentation XMM-SOC-CAL-TN-0083, on 
http : //xmm. vilspa. esa. es/docs/documents/). 

The RGS data showed an absorption feature at ~ 1.4 
keV. In case we fit this absorption line with a GABS com- 
ponent, the best-fitting line energy is 1.42 ± 0.06 keV, line 
width a = 70 ± 10 eV, and EW of ~ 35 and ~ 62 eV for 
RGS1 and RGS2, respectively. The width and strength of 
this absorption feature suggest that it is most likely due to 
a mismatch in the spectrum curvature of the PN and RGS 
data. In the rest of the analysis we fit this absorption line 
with a GABS component. 

Further, the EPIC fit residuals above 6 keV showed ad- 
ditional features in the PN spectrum at high energy (see 
Fig. |4|). The slight increase at ~ 8.3 keV (and the conse- 
quent dip at higher energies) is probably due to the emission 
lines in the detector background (see Fig. lA7[) . Since the ab- 
sorption feature at ~ 8.7 keV is only apparent in the PN 
data, and not in the MOS2 or PCA data, it is most likely 
an effect of not correcting for the background emission lines 
in the PN detector. To improve the fit statistics, we fit the 
absorption line at 8.78 ± 0.07 keV, with width a = 70 ± 20 
eV and EW of 37 ± 10 eV. 

In contrast to all above mentioned absorption features, 
the dip at ~ 9.2 keV could be real since both PN and PCA 
spectra show this feature. Since the highly ionised Fe line 
should be accompanied by an associated edge, we fit this 
feature with an edge component, and found a best-fitting 
edge depth of r ma x = 0.20 ± 0.02 at an energy of O.OSlp'Jy 
keV, midway between the He and H-like edges of Fe at ~ 
8.83 and ~ 9.28 keV, respectively. For PCA data only, the 
best-fitting T max = 0.19 ± 0.04 at an energy of 9.14 ± 0.17 
keV, consistent with the edge of H-likc Fe. 
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3.3.2 Line and Continuum models 

For testing different line and continuum models, we si- 
multaneously fit PN in the 1.3-11 keV range (excluding 
the 1.75-2.25 keV range), RGS in the 0.4-2 keV range, 
PCA in the 7-20 keV range, and HEXTE in the 20-150 
keV range. From the available XMM-Newton /EPIC data 
we only included PN data since this is the best calibrated 
instrument in timing mode, but for a consistency check we 
also fit the MOS2 data separately (see Section \3. 3.311 . 

Different relativistic line profiles: 

The broad and skewed Fe emission line profile that is 
apparent in the spectra of XTE J1652— 453 (see Fig. is 
consistent with what is expected from relativistic effects, and 
we therefore fit it with the commonly used line profiles for a 
non-rotating fpiSKLlNE; lFabian et al.ll"l989r i and a maximally 
rotating black hole (Laor). Both line components have sim- 
ilar parameters: line energy, El, which depends on the ion- 
isation stage of Fe, inner disc radius, Rj n , outer disc radius, 
.R ou t, disc inclination, 6, and emissivity index, /3, which is 
the power-law dependence of the disc emissivity that scales 
as R~P . The normalisation of the line, JVl, gives the number 
of photons cm -2 s _1 in the line. We let all parameters free, 
except i?out which we fixed to 400 GM/c 2 (with G and c 
the usual physical constants, and M the mass of the black 
hole). We further constrained the line energy to the range 
between 6.4-6.97 keV, since this is the expected range for 
neutral to hydrogen-like Fe-Ka lines. 

The fit with the diskline component did not provide 
a good fit with xt ~ 1-5 for 362 d.o.f., with strong residu- 
als below 6 keV, and the inner radius pegged at 6 GM/c 2 . 
Since the diskline appears not to be the right model to 
describe the line, we do not further mention any results for 
this model. The Laor profile however, gives an acceptable 
fit with xi ~ 1-2 for 362 d.o.f. In Table[T]we show the best- 
fitting parameters for the model with Laor, with a diskbb 
and powerlaw as our continuum. In the rest of the pa- 
per we refer to the phabs(diskbb + powerlaw + Laor) 
model as Model 1. The continuum parameters are the power- 
law photon index, V, power-law normalisation, iVr (in pho- 
tons keV -1 cm" 2 s" 1 at 1 keV), equivalent hydrogen column 
density, JVh, temperature at the inner disc radius, fcTin, and 
the disc blackbody normalisation, iVdbb = [(Ri n /km)/(D/10 
kpc)] 2 cos#, where Rj n is the apparent inner disc radius, D 
the distance to the source, and 8 the angle of the disc. 

The broadband spectrum of XTE J1652-453 and the fit 
residuals for Model 1 are shown in the left panel of Fig. O 
with in the inset the residuals in the Fe line region when 
fitting the underlying continuum spectrum with Model 1 
(except the line component) and excluding the Fe line band, 
i.e., the 4-7 keV range. In the right panel of Fig. [5]we show 
the unfolded spectrum and the model components of Model 
1. 

Since the Laor profile assumes a maximally spinning 
black hole, we also tried a relativistic line profile that al- 
lows to fit for the black hole spin (kyrline, up-dated ver- 
sion, private commu nication with M. Dovciak, but see also, 
iDovciak et aUl2004h ■ Besides the black hole angular momen- 
tum, a/M, the kyrline parameters are 6, Ri n , Rout, -El, 
and JVl, which are the same parameters (and same units) 
as for LAOR. We fixed -R ou t to 400 GM/c 2 and left all other 



parameters free. Further, within kyrline it is possible to 
set the directionality law for the emission (e.g., isotropic or 
Limb darkening/brightening), and to have a broken power- 
law radial emissivity profile. We assumed a disc with a single 
emissivity index, where j3 was free to vary. We further as- 
sumed a transfer function including Laor's limb darkening, 
and we only included emission from outside the marginally 
stable orbit. We fit kyrline using the same continuum as 
before, which is our Model 2: phabs(diskbb + powerlaw 
+ kyrline) . The best-fitting parameters for this model are 
also given in Table [T] 

From Table Q] (columns 1 and 2) and Fig [S] we can see 
that the spectrum of XTE J1652— 453 is well described 
by three components. The spectrum shows a significant 
contribution of thermal emission coming from a disc with a 
temperature of ~ 0.6 keV, where the disc fraction is ~ 75 
per cent. The spectrum has significant hard emission up to 
at least 40 keV, where the hard tail is described by a power 
law with a photon index of ~ 2.2. We further find that that 
the column density of interstellar matter in the direction of 
XTE J1652-453 is ~ 7 x 10 22 cm" 2 , which absorbs most of 
the X-rays coming from the disc. From Fig. [5] it is clear that 
there is a strong asymmetric Fe emission line with a red 
wing extending down to about 4 keV. From fits with the 
relativistic line profiles we find a line that is more than 15<r 
significant (based on the la error of the line flux), and a 
best-fitting line energy of 6.97 keV. Despite the fact that the 
kyrline model is more complex (having more parameters 
to fit for) and has a better grid resolution (more steps over 
which the line is integrated) than the Laor component, the 
inferred line parameters for both profiles are still consistent. 
From both relativistic line profiles we find that the disc is 
observed at a low inclination, < 20°, and has an inner disc 
radius of ~ 3.5 GM/c 2 . From the kyrline profile we find 
that the Fe emission line is indicative of a fast spinning 
black hole, although it is also consistent with a moderately 
spinning black hole, with a spin parameter of ~ 0.45. 

Different continuum models: 

As shown in previous work (see e.g., iHiemstra et al.l 

2009), the choice of continuum model can affect the param- 
eter values of the line and other continuum components, 
affecting any physical interpretation. To investigate the 
effect of the continuum on the inferred line profile, we 
tested different continuum components. Given that the 
best-fitting line parameters are consistent between the fits 
with a Laor or a kyrline component, we used the Laor 
profile to model the Fe emission line since it gives a slightly 
better fit, and it is computational much faster than kyrline. 

/: CUTOFFPL 

We first tried to fit the hard emission with a CUTOFFPL 
component to see whether there is an exponential roll-off 
at high energy. However, we did not find a roll-off in the 
hard power-law tail within the energy coverage of our data. 
Therefore, using a CUTOFFPL is effectively the same as 
using a powerlaw component. 

II: CompTT 

Next (Model 3), we replaced the powerlaw component in 
Model 1 with a component that describes the hard emis- 
sion as thermal Comptonisation of soft photons in a hot 
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Figure 5. Left: XMM-Newton / RXTE spectra of XTE J1652-453: RGS (0.4-2 keV, blue), PN (1.3-11 keV, black), PCA (7-20 keV, 
red), and HEXTE (20—150 keV, green) with, in the bottom panel, the best-fitting residuals of Model 1. In the sub-panel in the upper-right 
corner we show the residuals in the Fe line region when fitting the underlying continuum spectrum with Model 1 (without Laor) and 
excluding the Fe line band, i.e., the 4-7 keV range. Right: Unfolded XMM-Newton/ RXTE spectra of XTE J1652-453 and the model 
components (Model 1). For plotting purposes we binned the HEXTE data by a factor or 4. 



plasma (CompTT; iTitarchukl ll994T ). where we assumed a 
disc geometry for the Comptonising plasma. The parameters 
of CompTT are the input soft photon (Wien) temperature 
To, the plasma temperature kT, the plasma optical depth 
r, and the normalisation A^comp- The best-fitting results for 
Model 3 are given in Table [1] 

We find that the spectrum of XTE J1652-453 is 
well described by a thermal component and a Compton- 
ising component, where the Comptonising plasma is hot 
(kT ~ 170 keV) and optically thin (r <C 1). Because of the 
low-energy cutoff at the seed photon energy, the unabsorbed 
0.5-10 keV flux of the CompTT component is ~ 2-3 times 
lower than the flux of the powerlaw component in Models 
1 and 2. Due to this cutoff, the hard power-law tail does 
not contribute to very low energies and causes the slightly 
lower value for JVh, and thus a slightly higher contribution 
of the disk component. The other parameters are consistent 
with the values obtained in Models 1 and 2. 

Ill: CompPS 

Subsequently, we tried a component that describes a hy- 
brid thermal/non-thermal C omptonising plasma (CompPS; 
iPoutanen fe Svenssonl Il996l ). This model component com- 
putes the Comptonisation spectra taking into account the 
geometry and optical depth of the corona, the spectral dis- 
tribution of soft seed photons, and the parameters of the 
electron distribution. The resulting spectrum includes the 
incident thermal emission, the Comptonized emission, and 
the Compton reflection. For the source of soft seed photons 
(the incident thermal emission) we assumed a multicolor 
disk blackbody, and we therefore did not include a sepa- 
rate diskbb component to the model. Since we only wanted 
to test how a hybrid plasma affects the line and continuum 
parameters, we did not include the Compton reflection in 
CompPS, but we still added a Laor component to fit the Fe 
line. For a spherical geometry of the corona, the best-fitting 
results for this model (Model 4) are given in Table[T] The free 



parameters are the electron temperature of the corona, kT e , 
seed photon temperature, kTdbb, optical depth of the corona, 
r, and normalization, Ncps, which has the same units as 
the diskbb normalization. We further fixed the power-law 
index of the energy distribution of the non-thermal electrons 
to 2, we set the minimum and maximum lorentz factors of 
the non-thermal electrons to 2 and 3000, respectively, and 
set the abundances to solar. 

We find that Model 4 fits the data well with xt ~ 1.2 
for 362 d.o.f. The hybrid Comptonising plasma has a 
best-fitting temperature of ~ 75 keV and is optically thin 
(r ~ 0.15). The hybrid plasma is cooler than the fully 
thermal plasma in Model 3, while the optical depth is a 
factor 2 larger. We further find that the parameters of 
the disk and line are hardly affected by the choice of the 
Comptonising component (non-thermal, Model 1, thermal, 
Model 3, or hybrid, Model 4), with values being consistent 
within errors. However, the equivalent width is ~ 10-20 per 
cent lower in Model 4 than in the other 2 models. 

IV: REFLIONX 

Since the line energy of the Fe emission line suggests a high 
ionisation stag e of Fe, we tried an ionised reflection model 
(reflionx; I Ross fe Fabiani r2005). This model describes the 
reflection off an ionised and optically-thick atmosphere of 
constant density (such as the surface of an accretion disc) 
that is illuminated by radiation with a power-law spec- 
trum. Since reflionx only provides the reflected emission, 
we included a powerlaw component to account for the 
direct emission. We coupled the photon index of the illu- 
minating power-law component in reflionx to the pho- 
ton index of the powerlaw. reflionx includes ionisation 
states and transitions from the most important ions, like O 
III- VIII, Fe VI-XXVI, and several others. To account for 
relativistic effects in the disc, we used a relativistic con- 
volut ion component (Kerrconv; iBrenneman fe Reynolds! 
2006). All together, Model 5 was phabs (diskbb + power- 
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law + Kerrconv*reflionx) model. For this continuum 
model the absorption component at ~ 8.7 keV was not sig- 
nificant anymore and we therefore took it out. The absorp- 
tion component at ~ 1.4 keV and the edge at ~ 9.2 keV were 
still included. The free parameters in the reflionx model 
are the ionisation parameter, £ (where £ = 4nF/n, with F 
the total illuminating flux and n the hydrogen number den- 
sity), and the normalisation, N re i. We initially also let the Fe 
abundance to be free, finding a best-fitting value of 4 times 
solar. This is unlikely for a low-mass X-ray binary, and we 
therefore fixed the Fe abundance to solar. For Kerrconv, 
we assumed a single emissivity index, which was free to vary. 
The black hole spin and inclination angle of the disc were 
also free parameters. We further fixed the inner radius of 
the disc to be equal to the radius of the marginally stable 
orbit, Bms, and the outer radius to 400 R mB . The best-fitting 
parameters for Model 5 are also given in Table [T] 

We find that Model 5 fits the overall continuum and 
Fe emission line quite well, except that a clear absorption 
feature remained in the residuals at ~ 7.1 keV, causing 
a slightly higher xt compared with Models 1-3. Although 
this absorption feature is not obviously apparent in the raw 
data, we checked whether it can be either a K-edge of neu- 
tral Fe (7.11 keV), a blue-shifted absorption line of highly 
ionised Fe, or a mismatch in the model used to describe 
the data. We tried to fit the residuals at ~ 7.1 keV with 
an EDGE component, but since the feature is narrow, an 
edge does not provide a good fit. When we fit the residu- 
als with a GABS component, the best- fitting absorption line 
has a central energy of 7.2 keV and a width of 0.1 keV. 
For the case of a blue-shifted Fe-XXVI Ly-a line, the blue 
shift corresponds to a velocity of ~ 9500 km s , which 
is significantly higher than the blu e shift of ~ 300- 1600 
km s" 1 observed in GRO J1655-40 (jMiller et al.ll2006h . For 
other Fe lines the blue shift would be even larger. Since re- 
flionx self-consistently incorporates Compton scattering, 
the Fe feature in this model has an extended blue wing up 
to ~ 8 keV, whereas the observed Fe emission line in the 
spectrum of XTE J1652— 453 shows a sharp drop at ~ 7 
keV. We therefore conclude that the negative residuals at 
~ 7 .1 keV are likely caused by a deficiency in the model 
(see iDone fe Gierlihskil I2006T ) . 

From Model 5 we find a high ionisation parameter, 
£ ~ 2900 erg cm s~ , indicating that the emission line is 
dominated by highly ionised Fe, which is consistent with 
the line energy of 6.97 keV that we inferred for the relativis- 
tic line profiles in Models 1-3. We also find a moderately 
spinning black hole with a spin parameters of ~ 0.5, consis- 
tent with the lower confidence value inferred from Model 2. 
Further, the best-fitting inclination is consistent with that 
inferred from Laor, while the best-fitting emissivity index is 
slightly lower than we found from Models 1-3. The POWER- 
LAW component in Model 5 only represents the direct emis- 
sion, and therefore its normalisation is ~ 2 times lower than 
for the other three models. All other parameters of Model 
5 are consistent with those inferred for the other models. 
The reflection fraction for Model 5, as estimated from the 
(0.5-150 keV) flux ratio between the reflected and incident 
emission, is ~ 90 per cent. 



3.3.3 Consistency check: MOS2 data 

For consistency, we verified our results obtained from fits to 
the RGS/PN / RXTE data (in this section simply referred to 
as PN data) with fits to the RGS/MOS2/.RA' T£ data (and 
hence referred to as MOS2 data). For comparison we only fit 
Model 1 to the MOS2 data, assuming the same abundances, 
cross sections, and energy ranges that we used to fit the 
PN data (see Section \3. 3 \ . We ignored the MOS2 channels 
between 1.75-2.25 keV for similar reasons as for PN, added a 
Gaussian absorption component to fit the ~1.4 keV feature 
seen in the RGS spectra, and fit the ~ 9.2 keV feature in 
MOS2 with an edge. 

Model 1 (including all above absorption features) fit 
the MOS2 data well with X l = 1-25 (458.0/366). The un- 
absorbed 0.5-10 keV flux is 2.3 x 10" 9 erg cm" 2 s" 1 , and 
the best-fitting parameters for the continuum components 
are N H = (6.77 ± 0.05) x 10 22 cm" 2 , kT in = 0.60 ± 0.01 
keV, iVdbb = 859^4?, T = 2.17 ± 0.05, and Nr = 0.13 ± 0.02 
photons keV" 1 cm" 2 s" 1 at 1 keV. For the Fe line the best- 
fitting values are E L = 6.97t° . 51 keV, R in = 2.8 ± 0.5 
GM/c 2 , /3 = 3.7 ± 0.4, 6 = 25lg 8 degrees, and N L = 
1.9±0.5xl0" 3 photons cm 2 s 1 in the line. The Fe emission 
line has an EW of 514~^],® eV consistent with that obtained 
from PN data. The edge at ~ 9.2 and the absorption feature 
at ~ 1.4 keV have parameters which are consistent with the 
values obtained from the PN data. 

We find that most of the continuum and line param- 
eters obtained from the MOS2 data are consistent (within 
errors) with the values obtained from the PN data. The only 
parameters that were not consistent are the diskbb normal- 
isation and the inclination. However, the inclination inferred 
from the MOS2 data is consistent (within errors) with the 
inclination deduced from the kyrline profile that we fit to 
the PN data. The diskbb normalisation obtained from fits 
to the MOS2 data is lower than that for PN, which is most 
likely due to the fact that we subtracted the background 
spectrum in the case of the MOS2 data, but we did not for 
the PN data. 



3.4 Other possible broadening mechanisms 

By using the Laor and kyrline components in our Models 
1-4 we are in fact assuming that the line is only broad- 
ened by relativistic effects. However, one should consider 
the possibility that other mechanisms contribute to the to- 
tal line width. Besides broadening due to thermal Compton- 
isation by the disk, as in reflionx (see Section l3.3.2p . an 
alternative mechanism could be Compton broadening due 
to el ectron scattering in the hot plasma above the accretion 
disc l|Czernv et al.lll99ll ). The width of a Compton scattered 
line depends on the temperature and the optical depth of 
the Comptonising plasma. For a single Compton scatter- 
ing o-e/E = [(2fcT c )/(ra c c 2 )] 1/2 , with <t e the line width, 
E the central line energy, kT e the electron temperature 
of the plasma, and m e a nd c the usual physical constants 
l|Pozdnvakov et al.l ll983V From fits including a CompTT 
component to describe the high-energy emission (Model 3), 
we find r ~ 0.06 and kT a ~ 170 keV (see Table [TJ, which 
means that the Comptonising plasma would only intercept a 
fraction of 0.06 of the line photons, and scatter these into a 
broader line of width ~ 5.5 keV. A fraction of 0.94 of the line 
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Table 1. Different line and continuum models fit to XMM-Newton / RXTE data of XTE J1652-453 



Model component 


Parameter 


Model 1 


Model 2 


Model 3 


Model 4 


Model 5 






6.73+0.01 


6.75+0.02 


5554-0.01 


6.71+0.02 


6.73+0.01 


DISKBB 


fcT in (kcV) 
^dbb 

Flux (10~ 9 ) 


0.59±0.01 

loisj: 3 ^ 

2.09±0.07 


0.59±0.01 
1051±40 
2.12±0.08 


0.60±0.01 
1090+^5 
, 07+O.2O 
"'-0.08 




0.59+0.01 

998+10 
2.08+0.02 


POWERLAW 


r 

Flux (10~ 9 ) 


2.23+0.03 
0.18+0.01 
0.74±0.06 


2.22±0.03 
0.18±0.01 
0.73±0.05 






2.16+0.02 
0.09+0.01 
0.39+0.04 


CompTT 


T (kcV) 
kT (kcV) 

T 

JVComp (lO" 4 ) 

Flux (10~ 9 ) 






n aa+ - 12 
°- 88 -0.03 

170+ 24 

— 51 
0.06±0.01 
3.25^ 
,7+0.77 


- 




CompPS 


kT e (keV) 
kT dbb (kcV) 

T 

NcPS 

T7»1 / 1 M — Q \ 

t lux (10 ) 


: 


: 




75.3+1.3 
0.58+0.01 
0.15+0.01 

1314+37 
-7O + 0.12 
2 ' 72 -0.03 


: 


Laor 


B L (kcV) 
fl in (GM/c 2 ) 

e (dog) 

JV L (lO" 3 ) 
EW (oV) 


6 97+° 
°- -0.02 

3.4±0.1 

3-«tS:S 

1.6±0.1 


— 


6 97+° 

3.4+0.1 
, c ^+0.27 

9.4±3.9 
1.6+0.1 


6 97+° 
°' M '-0.02 

3.3 + 0.1 

401 + 0.18 

- 0+3.3 

1.3 + 0.1 

366+" 

— 15 


— 


KYRLINE 


a/M (GM/c) 
B L (kcV) 

B 

R ia (GM/c 2 ) 
6 (dog) 
N-L (10 -3 ) 
EW (oV) 


_ 

- 


0.96+°°] 

— 0.51 

6 97+° 
°- a '_0.04 

3.5±0.1 
3 39+ ' 31 
16.6±3.0 
1.8±0.2 
4 70 +28 

— 41 


_ 
- 


_ 
- 


_ 

- 


REFLIONX 


£ (erg cm s~ 1 ) 

Jv rof (10- 7 ) 










2872+* . 7 
4 g+0.04 
^-0.23 


Kerrconv 


a/M 

/3 

9 (dog) 










0.45+0.02 
2.7+0.1 
8.8+0.1 




xl (x 2 /v) 

Total Flux (10~ 9 ) 


1.22 (422.2/362) 
2.84±0.04 


1.25 (451.1/361) 
2.87±0.03 


1.21 (437.2/360) 

2 65+ - 76 
^■"^-o.og 


1.23 (445.2/362) 
2.73+0.02 


1.26 (460.2/365) 
2.86+0.11 



NOTES.- Results are based on spectral fits to RGS/PN/PCA/HEXTE data in the 0.4-150 keV range. Columns 1-2 give the results of 
different relativistic line models using the continuum model PHABS( DISKBB + POWERLAW) with line profiles Laor (Model 1) and KYRLINE (Mo- 
del 2). Columns 3—5 give the results of different continuum models. Models 3 and 4, respectively, include a thermal Comptonising plasma 
(CompTT) and a hybrid plasma (CompPS) to model the hard emission, and the Laor component is used to fit the Fc line. For Model 4 the 
Comptonising spectrum sclf-consistently includes the incident thermal emission, so we did not include a separate DISKBB component as we 
did for Model 3. Model 5 is the continuum model phabs(diskbb + POWERLAW + Kerrconv*reflionx), where REFLIONX describes the reflected 
emission, POWERLAW models the direct emission, and Kerrconv smears for relativistic effects. For details about parameters, see text. Cross- 
calibration constants for the different instruments with respect to PN are: 1.09 (PCA), 0.84 (RGS1 & HEXTE), and 0.79 (RGS2). Fluxes 
are unabsorbed and given for the 0.5—10 keV range in units of erg cm 2 s 1 . 



is unaffected, and comes out at its intrinsic width. Thus the 
line would have a predominantly narrow core, but with 6 per 
cent of the photons scattered to such extent that it would 
be difficult (or even impossible) to distinguish from the con- 
tinuum. However, as a matter of fact the problem here is 
that in all probability the Comptonisation is not thermal 
but instead is non-thermal. Therefore the obtained optical 
depth and temperature are not real. As we found from the 
fits of a hybrid Comptonisation model (Model 4) , the plasma 
can also be cooler (kT e ~ 75 keV), but we still found a low 
optical depth (r ~ 0.15). In this case ~ 15 per cent of the 
line photons will be scattered into a broader line of width 



~ 3.5 keV, which is still too broad and too weak to be dis- 
tinguished from the continuum. 

iKallman fc White! (QUI) showed that the apparent 
widths of Fe lines may be due to a blend of narrow line 
components arising from different ionisation states, with 
each of these narrow components in turn affected by Comp- 
ton broadening. The combination of blending and Compton 
broadening produces a line profile with an equivalent width, 
line width and line energy that depend on several factors, 
besides the ionisation parameter of the emitting gas. In this 
scenario, the total line width is not only set by relativis- 
tic effects, but also by the combined effects of gas motions 
(rotation), Compton scattering, and blending of emission of 
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several ionisation states, where the centroid line energy is 
a function of the degree of ionisation of Fe, and hence the 
ionisa tion parameter of the emitting gas (jKallman fc White! 
1989). To estimate the effect of line blending, we fit the Fe 
line with a multi-component line model. Assuming the con- 
tinuum of Model 1, we added three Laor components with 
the line energies fixed at 6.4, 6.7 and 6.97 keV, for neutral, 
He- like, and H-like Fe, respectively. We fixed the outer radius 
of the disc to 400 GM/c 2 , and let the rest of the parameters 
free to vary. Except for the line normalisations, we coupled 
the line parameters of each of the three Laor profiles. We 
find that the multi-line model fits the observed Fe emission 
line well, xt = 1-17 for 361 d.o.f. The best-fitting parame- 
ters are R in = 2.201^° GM/c 2 , /3 = 4.0±0.1, 9 = 40.7±1.5 
degrees, and iV L = 2.0 ± 0.2 x 10~ 3 , and 0.6 ± 0.28 x 10~ 3 
photons cm -2 s _1 for the 6.4 and 6.97 keV lines, respec- 
tively. The 6.7 keV line was not required by the fit, with UL 
to the EW of ~ 51 eV. The equivalent widths of the neutral 
and H-like Fe were 389tf 9 and 127±41 eV, respectively. The 
best-fitting continuum parameters were consistent with the 
values reported for a single- line model (Model 1, Table [T]). 
From this we initially concluded that line blending can par- 
tially explain the observed line profile, although relativis- 
tic effects are still important. The resulting disc inclination 
was significantly higher for the multi-line model than for 
the single-line model. In contrast to the single Laor profile, 
the multi-line model was dominated by the neutral-Fe line, 
with an equivalent width ~ 3 times higher than the one of 
the H-like Fe line. Based on the high ionisation parameter, 
£ ~ 2900 erg cm s _1 , obtained from the fits with reflionx 
(see Table [T]) we do not expect a significant fraction of neu- 
tral Fe in the disc, which casts some doubt on this approach. 
Moreover, given the observed high disc temperature we do 
not expect much neutral Fe, irrespective of photo-ionisation. 

A third broadening mechanism that we considered ex- 
plains the red-skewed emission line by down-scattering of 
line photons in an optically thick and par tly ionised wind 
from the BH (see iTitarchuk et al.l 120091 . and references 
therein). In this model the effect of photo-absorption of the 
line photons at higher energies (above 7-8 keV, depending 
on the ionisation stage of the wind) suppresses the blue wing 
of the Fe emission line, and therefore leads only to a red- 
shifted profile, where the redshift of the line photons is to 
first order a function of v/c, with v the outflow velocity. 
Contrary to the optically thick Comptonising wind assumed 
in this model, we find an optically thin corona for the case of 
XTE J1652— 453, which makes this outflow model unlikely 
to explain the observed broad and red-skewed Fe line. 

3.5 Timing analysis 

The power spectrum of the RXTE observation that was 
taken simultaneously with the one of XMM-Newton,can 
be reasonably fit (\ 2 = 35.9 for 25 d.o.f.) with three 
Lorentzians at 2.2 ± 0.4 Hz, 10.0 ± 1.7 Hz, and 131 ± 20 
Hz, respectively. (These are the frequencies at which the 
Lore ntzians peak in the frequency x power representation of 
e.g.. iBelloni et aLll2002T ) The single-trial significance of each 
Lorentzian (defined as the power of the Lorentzian divided 
by its negative error) was 4.9cr, 2.8a, and 3.4a, respec- 
tively. The three Lorentzian had, respectively, quality fac- 
tors, Q = fo/FWHM (where v Q and FWHM are the cen- 




v (Hz) 

Figure 6. Power density spectrum of XTE J1652— 453 using 
the RXTE /PCA data from the observation taken simultane- 
ously with that of XMM-Newton (MJD 55065, RXTE ObsID 
94432-01-04-00). The plot shows the frequency x Fourier Power in 
units of squared rms fractional amplitude versus frequency. The 
black solid line is the best-fitting model to the data using three 
Lorentzians. The three Lorentzians are shown by the red dashed, 
green dot-dashed, and blue dot-dot-dashed lines, respectively. 



tral frequency and the full-width at half maximum of the 
Lorentzian) of 0.5 ± 0.2, 1.0 ± 0.7, and 1.5 ± 0.9, and total 
fractional rms amplitudes of 8.9T0.9, 8.3T1.6, and 13.8T2.2 
per cent. In Fig. [6] we show the power density spectrum of 
this observation, together with the best-fitting model. The 
black solid line is the total model, while the dashed, dot- 
ted and dashed-dotted lines in colour show the individual 
components. 

If the marginally significant Lorentzian at ~ 130 Hz is 
real, we can use it to put an upper limit on the black-hole 
mass, as any frequency should be less than Keplerian. For a 
Keplerian frequency of v ~ 130 Hz at the inner edge of the 
accretion disc of R in ~ 4 GM/c 2 , v = c 3 /[2nGM(Rf^ 2 + a)] 
|Bardeen et alJI 19721 ). with a the spin parameter and M the 
mass of the black hole, we find an upper limit on the mass 
of ~ 30Mq . Given that most known BH binaries have mass 
around ~ WMq , then this shows that these fluctuations are 
indeed produced close to the black hole. 



4 DISCUSSION 

The XMM-Newton/ RXTE energy spectrum of XTE 
J1652— 453 shows a broad and strong Fe emission line with 
an equivalent width of ~ 450 eV, which is among the 
strongest ever se e n from a black-ho le candidate (see e.g., 
iRossi et al.l 120051 ; iMiller et all 120041 '). The line is consis- 
tent with being produced by reflection off the accretion 
disc, broadened by relativistic effects. The power spectrum 
from RXTE data can be fitted with two Lorentzians be- 
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low 10 Hz, and shows an additional variability component 
above 50 Hz. The simultaneous XMM-Newton / RXTE ob- 
servations were performed when XTE J1652— 453 was in 
the hard-intermediate state ( HIMS; iMendez fc van der Klisl 
1 19971 ; iHoman fc Bellonil fcoOS ) , as the outburst decayed from 
the disc dominated high/soft state towards the low/hard 
state. Our high-resolution, broadband observations substan- 
tially increase the number of good datasets from this state 
(|Miller et al.l |2004 the other is GX 339-4). The strong 
line as well as the additional component at high frequen- 
cies are both char acteristic of the hard-intermediate state. 
Miller et al.l (|2004T ) show a Chandra grating spectrum of the 
black-hole candidate GX 339—4, with a very similar contin- 
uum spectrum as XTE J1652— 453, where the line has an 
equiva lent width of ~ 600 eV. iKlein-Wolt fc van der Klisl 
(2008) show that there is an additional noise compo- 
nent at high frequencies in the HIMS, although it is 
somewhat weaker than the component seen in our data. 
Power spectra taken in the HIMS also typically show 
the strongest, most coherent, low-frequency QPO (e.g., 
IKlein-Wolt fc van der Khsil200Sft . The lack of any such co- 
herent low-frequency QPO in our data may imply that we 
view the system at a rather low inclination angle. This is 
consistent with the inclination inferred from fitting the Fe 
line profile with relativistic models. 

Understanding the intermediate state is important as 
it marks the onset of the spectral transition, and the ap- 
pearance of the radio jet. While there is still debate on the 
nature of the transition, most models currently explain it 
by t he inner disc evaporating into a h ot coronal flow (see 
e.g.. iLiu etal]|2006l ; iMever et al.l 120001 ) . The intermediate 
state is therefore the one in which the disc is only slightly 
recessed from the last stable orbit, so there is strong disc 
emission, but there is also strong hard X-ray flux from the 
re-emerging hot flow. This geometry also predicts that there 
should be strong reflection features, and that these should be 
strongly smeared by relativistic effects. As the outburst de- 
clines further, the models predict that the inner disc quickly 
recedes, being replaced by the hot flow. With the receding 
disc, the solid angle subtended by the optically thick disc 
decreases, the amount of reflection decreases, and the larger 
radius means that the line is less strongly smeared. These 
models then suggest that the strongest and broadest Fe lines 
should be seen in the intermediate states, where the disc is 
closest to the last stable orbit, but co-exists with enough 
of the hard X-ray emitting corona to produce the required 
illumination of the disc. 

From low-resolution (RXTE) spectral monitoring it is 
indeed seen that for the high/soft states, as there is very 
little hard X-ray emission illuminating the disc, the Fe lines 
are generally weak, while they appear strongly at the in- 
termediate states, as coronal X-ray emission re-appear. The 
line wi dth and streng th then decrease as the outburst de- 
clines |Gilfanovl [20101. However, such low-resolution data 
cannot distinguish changes in the reflection broadening from 
changes in the ionisation state, and since the disc tempera- 
ture and/or strength of the irradiation change, the ionisation 
state should also change. Even if there is no change in ge- 
ometry, an increase in disc temperature gives an increase in 
ionisation state, and so the strength of the conseque nt reflec- 
tion features increases (Ross fc FabiarJ[l993l , l2007f) . as well 
as the line width increases due to increased Compton broad- 



ening (e.g.. lYoung et alj|2005ft . Moderate to high resolution 
data show that the Fe lines are ionised in the intermedi- 
ate state (e.g., I Miller et al. 1 12004 and our data), and are 
predominantly neutral in the hardest low /hard state (e.g., 
iTomsick et~aT]|2009l ; [Hiemstra et alj|2009h . 

The high-resolution, broadband spectrum of XTE 
J1652— 453 allows us to separately constrain ionisation, solid 
angle and relativistic smearing. The ionisation is high with 
log£ ~ 3.5, and the very broad Fe line requires strong 
relativistic effects. The fits using a self-consistent reflec- 
tion model give Ri n — 4 GM/c 2 corresponding to the 
last stable orbit for a black hole of spin ~ 0.5. This 
is a lower limit since we expect the disc to be slightly 
recessed. None the less, the RXTE monitoring of XTE 
J1652— 453 shows that the direct disc emission has sim- 
ilar temperature and normalisati on as in the high/soft 
state earlier in the outburst ( Markwardt et al. 2009al ; 
iMarkwardt et al. 200 9b). supporting the idea that the disc 
is very close to the last stable orbit in our data. 

However, the fits with the self-consistent reflection 
model show clear negative residuals at ~ 7.1 keV. This 
could correspond to a blue-shifted absorption line, such as 
the ones that are oft en seen in high inclination systems (e.g., 
iKubota et al.ll2007T ) . This absorption feature is not obviously 
apparent in the raw data, and it is only required when the 
smeared reflionx model is used since the intrinsic Comp- 
ton broadening in this model predicts no sharp drop on 
the blue wing of the line (|Done fc Gier lihski 200(1). The re- 
flionx model underestimates the Compton broadening as 
it calculates the temperatur e of the disc from i rradiation 
only, for kT of 0.1-0.3 keV l|Ross fc Fabianl 1 1993ft . This is 
obviously lower than the observed inner disc temperature 
of 0.6 keV. The Compton broadened line width should be 
a oc E lincy /(2kT e )/(m e c 2 ) ~ 0.32 keV for He-like Fe. This 
is negligible compared to the 1.7-2 keV broadening implied 
by -Rin ~ 4 GM/c 2 , and should not significantly change the 
derived inner radius. The Compton broadening smooths out 
the sharp drop from the blue wing of the line even more, and 
increases the mismatch with the data at ~ 7.1 keV. Since 
the model does not describe the data in detail, we should be 
cautious in interpreting the reflection parameters in detail. 



5 SUMMARY 

We observed the black-hole candidate and X-ray transient 
XTE J1652-453 simultaneously with XMM-Newton and 
RXTE. During this observation the source w as in the hard- 
intermediate state (|Homan fc Bellonil |2005ft . We found a 
strong and asymmetric Fe emission line in the X-ray spec- 
trum of XTE J1652-453, with an equivalent width of ~ 450 
eV, one of the strongest lines ever seen from a black-hole 
candidate. From fits to only XMM-Newton /PN spectrum, 
the equivalent width is slightly lower, ~ 380 eV. 

The line is also broad, giving an inner disc radius of 
~ 4 GM/c 2 , from either a single Laor line profile or from a 
self-consistent ionised reflection model. This implies a spin 
parameter of ~ 0.5, which is a lower limit as the disc may 
be slightly recessed in this state. In any case, neither model 
fits the data well in detail. This is not at all surprising for 
the single Laor line as the line must be accompanied by 
a reflected continuum. Be that as it may, the self-consistent 
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ionised reflection model used here incorporates all the known 
effects such as multiple lines and edges from the different 
ionisation states of Fe (and other elements ), together with 
Compton broadening (|Ross fe Fabianl 12005). such that the 
line is intrinsically broad. The characteristic skewed profile 
expected from a single narrow line is smeared out, so there 
is no sharp drop in the model at the blue side of the line. 
However, the data clearly prefer such a feature, indicating 
that even the best current reflection models may not be 
able to describe the detailed shape of the Fe line seen in the 
spectrum of this state. 

The Fourier power spectrum of the RXTE observa- 
tion simultaneous with that of XMM-Newton shows addi- 
tional variability above 50 Hz. No coherent QPOs at low 
frequency are accompanying this, which may suggest that 
we see the source at a fairly low inclination angle. The 
3.4<r (single trial) component at ~ f30 Hz is too broad 
to be the high-frequency QPO, but in case it is, and as- 
suming that it reflects the Keplerian frequency of the in- 
ner disc radius of ~ 4 gravitational radii (inferred from the 
line profile), then the estimated upper limit on the black- 
hole mass is ~ 30 Mq. None the less, since the source 
is in the hard-intermediate state rather than in the soft- 
intermediate st ate where the high-frequency QP Os are oc- 
casionally seen l|Klein-Wolt fe van der Klis|[2008l ), this addi- 
tional high-frequency component may be a related feature 
at lower frequencies and lower coherence, and it might be 
connected to the re-emergence of the jet. 



6 ACKNOWLEDGEMENTS 

We thank the XMM-Newton team and the project scien- 
tist Norbert Schartel for performing the ToO observation. 
We also like to thank the RXTE team for rescheduling to 
obtain a contemporaneous observation. We are grateful to 
an anonymous referree for constructive comments, and we 
thank Juri Poutanen for a valuable discussion. PC acknowl- 
edges funding via a EU Marie Curie Intra-European Fellow- 
ship under contract no. 2009-237722 



REFERENCES 

Anders E., Grevesse N., 1989, Geochim. Cosmochim. Acta, 
53, 197 

Arnaud K. A., 1996, in G. H. Jacoby & J. Barnes ed., Astro- 
nomical Data Analysis Software and Systems V Vol. 101 
of Astronomical Society of the Pacific Conference Series, 
XSPEC: The First Ten Years, p. 17 

Bardeen J. M., Press W. H., Teukolsky S. A., 1972, ApJ, 
178, 347 

Belloni T., Homan J., Casella P., van der Klis M., Nespoli 
E., Lewin W. H. C, Miller J. M., Mendez M., 2005, A&A, 
440, 207 

Belloni T., Psaltis D., van der Klis M., 2002, ApJ, 572, 392 
Belloni T. M., 2010, in T. Belloni ed., Lecture Notes in 

Physics, Berlin Springer Verlag Vol. 794 of Lecture Notes 

in Physics, Berlin Springer Verlag, States and Transitions 

in Black Hole Binaries, p. 53 
Boirin L., Mendez M., Diaz Trigo M., Parmar A. N., Kaas- 

tra J. S., 2005, A&A, 436, 195 



Brenneman L. W., Reynolds C. S., 2006, ApJ, 652, 1028 
Cabanac C, Fender R. P., Dunn R. J. H., Kording E. C, 

2009, MNRAS, 396, 1415 
Calvelo D. E., Tzioumis T., Corbel S., Brocksopp C, 

Fender R. P., 2009, The Astronomer's Telegram, 2135 
Coriat M., Rodriguez J., Corbel S., Tomsick J. A., 2009, 

The Astronomer's Telegram, 2219 
Czerny B., Zbyszewska M., Raine D. J., 1991, in A. Treves, 

G. C. Perola, & L. Stella ed., Iron Line Diagnostics in X- 

ray Sources Vol. 385 of Lecture Notes in Physics, Berlin 

Springer Verlag, Compton Broadening of the Iron Line in 

NGC 3227. p. 226 
Done C, Diaz Trigo M., 2009, ArXiv e-prints 
Done C, Gierlihski M., 2006, MNRAS, 367, 659 
Done C, Gierlinski M., Kubota A., 2007, A&A Rev., 15, 1 
Dovciak M., Karas V., Yaqoob T., 2004, ApJS, 153, 205 
Fabian A. C, Rees M. J., Stella L., White N. E., 1989, 

MNRAS, 238, 729 
Fabian A. C, Zoghbi A., Ross R. R., Uttley P., Gallo L. C, 

Brandt W. N., Blustin A. J., Boiler T., Caballero-Garcia 

M. D., Larsson J., Miller J. M., Miniutti G., Ponti C, 

Reis R. C, Reynolds C. S., Tanaka Y., Young A. J., 2009, 

Nature, 459, 540 
Gilfanov M., 2010, in T. Belloni ed., Lecture Notes in 

Physics, Berlin Springer Verlag Vol. 794 of Lecture Notes 

in Physics, Berlin Springer Verlag, X-Ray Emission from 

Black-Hole Binaries, p. 17 
Hiemstra B., Soleri P., Mendez M., Belloni T., Mostafa R., 

Wijnands R., 2009, MNRAS, 394, 2080 
Homan J., Belloni T., 2005, Ap&SS, 300, 107 
Jahoda K, Markwardt C. B., Radeva Y., Rots A. H., Stark 

M. J., Swank J. H., Strohmayer T. E., Zhang W., 2006, 

ApJS, 163, 401 
Kallman T., White N. E., 1989, ApJ, 341, 955 
Klein- Wolt M., van der Klis M., 2008, ApJ, 675, 1407 
Kubota A., Dotani T., Cottam J., Kotani T., Done 

C, Ueda Y., Fabian A. C, Yasuda T., Takahashi H., 
Fukazawa Y., Yamaoka K, Makishima K, Yamada S., 
Kohmura T., Angelini L., 2007, PAS J, 59, 185 

Laor A., 1991, ApJ, 376, 90 

Lasota J., 2001, New Astronomy Review, 45, 449 
Liu B. F., Meyer F., Meyer-Hofmeister E., 2006, A&A, 454, 
L9 

Maccarone T. J., 2003, A&A, 409, 697 

Markwardt C. B., Beardmore A. P., Miller J., Swank J. H., 

2009, The Astronomer's Telegram, 2120 
Markwardt C. B., Swank J. H., 2009, The Astronomer's 

Telegram, 2108 
Markwardt C. B., Swank J. H., Krimm H. A., Pereira 

D. , Strohmayer T. E., 2009, The Astronomer's Telegram, 
2107 

Martocchia A., Matt G., Belloni T., Feroci M., Karas V., 
Ponti G., 2006, A&A, 448, 677 

McClintock J. E., Remillard R. A., 2006, Black hole bina- 
ries, pp 157-213 

Mendez M., van der Klis M., 1997, ApJ, 479, 926 

Meyer F., Liu B. F., Meyer-Hofmeister E., 2000, A&A, 361, 
175 

Miller J. M., 2007, ARA&A, 45, 441 

Miller J. M., Raymond J., Fabian A., Steeghs D., Homan J., 
Reynolds C, van der Klis M., Wijnands R., 2006, Nature, 
441, 953 



© 2010 RAS, MNRAS 000,rTlfT5l 



14 B. Hiemstra et al. 



Miller J. M., Raymond J., Fabian A. C, Homan J., Nowak 
M. A., Wijnands R., van der Klis M., Belloni T., Tomsick 
J. A., Smith D. M., Charles P. A., Lewin W. H. G., 2004, 
ApJ, 601, 450 

Mitsuda K., Inoue H., Koyama K., Makishima K., Mat- 
suoka M., Ogawara Y., Suzuki K., Tanaka Y., Shibazaki 
N., Hirano T., 1984, PASJ, 36, 741 
Poutanen J., Svensson R., 1996, ApJ, 470, 249 
Pozdnyakov L. A., Sobol I. M., Syunyaev R. A., 1983, As- 
trophysics and Space Physics Reviews, 2, 189 
Remillard R. A., McClintock J. E., 2006, ARA&A, 44, 49 
Reynolds C. S., Fabian A. C, 2008, ApJ, 675, 1048 
Reynolds M., Callanan P., Nagayama T., 2009, The As- 
tronomer's Telegram, 2125 
Ross R. R., Fabian A. C, 1993, MNRAS, 261, 74 
Ross R. R., Fabian A. C, 2005, MNRAS, 358, 211 
Ross R. R., Fabian A. C, 2007, MNRAS, 381, 1697 
Rossi S., Homan J., Miller J. M., Belloni T., 2005, MNRAS, 
360, 763 

Sala G., Greiner J., Ajello M., Primak N., 2008, A&A, 489, 
1239 

Sala G., Greiner J., Haberl F., Kendziorra E., Dennerl 
K., Freyberg M., Hasinger G., 2006, in A. Wilson ed., 
The X-ray Universe 2005 Vol. 604 of ESA Special Publi- 
cation, XMM-Newton Observations of the Microquasars 
GRO J1655-40 and GRS 1915+105. p. 291 

Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 

Thorne K. S., Price R. H., 1975, ApJ, 195, L101 

Titarchuk L., 1994, ApJ, 434, 570 

Titarchuk L., Laurent P., Shaposhnikov N., 2009, ApJ, 700, 
1831 

Tomsick J. A., Yamaoka K., Corbel S., Kaaret P., Kalemci 

E., Migliari S., 2009, ApJ, 707, L87 
Torres M. A. P., Steeghs D., Jonker P. G., McClintock J. E., 

Morrell N., Roth M., 2009, The Astronomer's Telegram, 

2190 

van der Klis M., 1989, in Ogelman H., van den Heuvel 
E. P. J., eds, Timing Neutron Stars Fourier techniques in 
X-ray timing, p. 27 

van Straaten S., van der Klis M., Mendez M., 2003, ApJ, 
596, 1155 

Verner D. A., Ferland G. J., Korista K. T., Yakovlev D. G., 
1996, ApJ, 465, 487 

Wilms J., Allen A., McCray R., 2000, ApJ, 542, 914 

Yamada S., Makishima K., Uehara Y., Nakazawa K., Taka- 
hashi H., Dotani T., Ueda Y., Ebisawa K., Kubota A., 
Gandhi P., 2009, ApJ, 707, L109 

Young A. J., Lee J. C, Fabian A. C, Reynolds C. S., Gib- 
son R. R., Canizares C. R., 2005, ApJ, 631, 733 

Zhang W., Jahoda K., Swank J. H., Morgan E. H., Giles 
A. B., 1995, ApJ, 449, 930 



APPENDIX A: CHOICE OF BACKGROUND 
AND ABUNDANCES 

Al Background in PN timing mode 

As can be seen in Fig. IA7I the PN background spectrum 
(red) extracted from the XTE J1655— 453 observation, is 
contaminated with source photons (see Section 12.2.11 for 
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Figure A7. Source and background spectra extracted from 
XMM-Newton data with PN in timing mode. The background 
and source spectra are all rebinned using the tool pharbn (see 
Section |2. 2, lj l. The black (top) spectrum is the source spectrum 
of XTE J1652-453, extracted from the RAWX columns 30 to 46; 
the red (middle) spectrum is the background spectrum of XTE 
J1652-453 using RAWX in [5:21]; and the blue (bottom) spec- 
trum is the background spectrum extracted from a 'blank-sky 
field' observation of 4U 1608—52 using the RAWX columns 10 to 
26. 



more details about the extraction criteria). Using this back- 
ground spectrum leads to an over-subtraction of the source 
spectrum and it may change the results of the spectral anal- 
ysis. To investigate whether the continuum and line parame- 
ters are affected by the background correction, we performed 
three fits to the source spectrum of XTE J1652— 453, using a 
different background correction in each fit. As a first back- 
ground correction (Model la) we used the background spec- 
trum extracted from the outer parts of the CCD of the same 
observation of XTE J1652-453 (RAWX in [5:21]). This is 
the background spectrum contaminated with source pho- 
tons. For the second background correction (Model lb) we 
extracted a spectrum from a so-called blank-sky field obser- 
vation that was taken with PN in timing mode, and as our 
last correction (Model lc) we did not subtract any back- 
ground. 

Since the background can be different in different parts 
of the sky, for a blank-sky field we chose an observation taken 
in a direction close to that of XTE J1652— 453; we therefore 
used the observation of the neutron-star low-mass X-ray bi- 
nary 4U 1608-52 (obsID 0074140201) taken just after an 
outburst, where the source had returned to quiescence. From 
all the blank-sky field observations taken in a direction close 
to XTE J1652-453, the one of 4U 1608-52 has the most 
similar Nn value compared with XTE J1652— 453. Since 
in the observation of 4U 1608—52 the source was still de- 
tectable, although it was weak, we extracted the background 
spectrum using the columns RAWX in [10:26], confirming 
that we do not include any columns contaminated with 
source photons. In Fig. IA7I we show the background spec- 
trum extracted from the 4U 1608—52 observation in blue. As 
can be seen from both background spectra in Fig. IA7I there 
are emission features at ~ 8 keV that are due to fluorescence 
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from the detectors and the structure surrounding them. The 
lines at ~ 8 keV are the fluorescence lines of Cu-Ka and 
Zn-Ka (see the XMM-Newton Users Handbook), and are 
time-dependent. Using a background spectrum from a dif- 
ferent observation than the source spectrum, as we do for 
4U 1608—52, may cause an incorrect background subtrac- 
tion at around ~ 7-8 keV and may change the parameters 
of a possible Fe line at ~ 7 keV. 

For the three spectral fits with different background cor- 
rections, we fit Model 1 to PN data only, covering the 1.3-10 
keV range, and excluding the 1.75-2.25 keV range. As de- 
scribed in Section T3. 31 we added an edge at ~ 9.2 keV. The 
results for the different background corrections are given in 
Table [AH 

We set the abundance s according to values reported in 
Anders fe Grevesse (|l989l ). and we used cross sections of 



Verner et al.l fl996J). For the purpose of comparing differ- 
ent background corrections, the assumed abundances and 
cross sections are not important. However, the best-fitting 
parameters and their physical interpretation can be affected 
when different abundance/cross-section tables are used (see 
following Section). For comparison, in Table [A"2l we also give 
the best-fitting re sults of Mode l lc, w here we instead used 
the abundances of I Wilms et"afl l|2000h (Model Id). 

We find that for Mod els la-c (i.e., with abundances 
of I Anders fe Grevessej [l989h all continuum parameters, ex- 
cept for JVh , are consistent between the three different back- 
ground corrections. The column density in Model la is 
slightly higher than in Models lb and lc. The diskbb nor- 
malisation is slightly lower in Model la than in Models lb 
and lc, but is consistent within errors. The higher JVh can be 
due to the over-subtraction of th e background, which also 
explains the lower JVdbb (see also ICabanac et al.l 120091 . for 
discussion about JVh versus JVdbb relation). We conclude that 
for the case of XTE J1652— 453 the choice of background cor- 
rection has no effect on the continuum and line parameters. 
However, the issue of a contaminated background would be 
relevant if the diskbb normalisation is used to determine 
the inner disc radius. 

The Xv given in Table lA"2l are relatively high, although 
the parameters obtained from fits to only PN data are con- 
sistent with the fits to XMM-Newton/ RXTE data, where 
the reduced \ 2 i s acceptable. The contribution to the high 
xt is most likely due to uncertainties in the detector cali- 
bration. An additional systematic error of 0.6 per cent to all 
PN channels is sufficient to get a xt ~ 1- 



fitting the continuum spectrum in a relatively narrow band 
(1.3-10 keV, only PN) or in a wide-coverage (0.4-150 keV) 
did not change the parameters significantly. However, fits 
to only PN data resulted in larger errors on the powerlaw 
parameters, and also resulted in a higher powerlaw nor- 
malisation (although still consistent within errors) than for 
the broadband fits. With the higher powerlaw normalisa- 
tion the equivalent width of the Fe emission line is clearly 
lower in the fits to PN data only than for the broadband fits. 
Nevertheless, for the case of XTE J1652— 453 we found that 
even for the different continuum parameters, the Fe line pa- 
rameters are not affected significantly. The best-fitting line 
parameters are consistent for the different background cor- 
rections, abundances, and energy coverage that we used. 

This paper has been typeset from a T[hX/ MpX file prepared 
by the author. 



A2 Abundances, cross sections, and energy 
coverage 

By investigating the effects of the different background cor- 
rections, we have used different abundances than in the 
rest of our analysis. By comparing the best-fitting results 
of Model Id with those of Models la-c in Table [SH we find 
that the equivalent hydrogen column density is ~ 30 per 
cent l ower if we used the abundances of lAnders fe Grevesse 
1 1989h than if we used the abundances of Iwilms et al. 
(2000). With the change in JVh, the diskbb normalisation 
is significantly affected, and with this the other continuum 
parameters are slightly changed. Further, the best-fitting pa- 
rameters of Model Id, although slightly different, are con- 
sistent with those for Model 1 (Table [TJ. In other words, 
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Table A2. Spectral fit results for different background corrections 



Model component 


Parameter 


Model la 


IVlOClCl ID 


Model lc 


IVloclcl Icl 


PHABS 


N H (10 22 ) 


4.75±0.03 


4.66±0.04 


4.64+0.03 


6.78+0.03 


DISKBB 


kT in (kcV) 
JVdbb 


0.56±0.01 
1299+.- 


0.56±0.01 
1428+gj 


0.56+0.01 
1390+.g 


0.58+0.01 
1060+. 2 * 


POWERLAW 


r 


2.52±0.09 
0.33+0.06 


2.52±0.09 
0.33+0.06 


2.47+0.05 
0.30+0.06 


2.34+0.07 
0.23+0.04 


Laor 


E L (kcV) 

P 

R in (GM/c 2 ) 

e (dog) 

AT L (IO- 3 ) 
EW (gV) 


6 97+° 
°- a -0.02 

3.3+0. 1 

u -0.05 

— o .0 
1.2±0.1 


6 97+° 

-0.03 

3.3+0.1 
401 + 0.29 
*- ui -0.46 

4 4+6-3 
-4.4 

1.2+0.3 

348 ^32 


6 97+° 
ua -0.04 

3.3+0.1 

1.2+0.1 

„.„ + 36 
3*3-29 


6 97+° 
°- a -0.02 

3.4+0.1 

401 + 0.19 

u -0.07 

5 1+ 2 - 9 
&1 -5.1 

1.3+0.1 

375+ 29 




xl (x 2 /-) 

Total flux (1U~ 9 ) 


1.58 (261.0/165) 
8.99 


1.49 (246.6/165) 
9.31 


1.45 (240.0/165) 
9.28 


1.50 (247.4/165) 
8.87 



NOTES.- Results are based on spectral fits to PN data of XTE J1652 — 453 in the 1.3-10 keV range. Background cor- 
rections are: background of XTE J1652 — 453 contaminated with source photons (Model la); background of 4U 1608 — 52 
(Model lb), and no background subtracted (Model lc 8z Id). The flux is unabsorbed and given for the 2—10 keV range, 
in units of erg cm - 2 s — 1 . Cross sections are set to^erner^t^L^ (^^96), and abundances to^^iders^^^Grevesse (^989^ 
for Models la— c, and to^ilm^^l^ ^2000^) for Model Id. Details about parameters are given in Section [373] 
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